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A unique feature of ultracold atoms is the separation of length scales, ro <g kf , where Zcf and ro 
are the Fermi momentum characterizing the average particle distance and the range of interaction 
between atoms respectively. For s-wave scattering, Shina Tan discovered that such diluteness leads 
to universal relations, all of which are governed by contact, among a wide range of thermodynamic 
quantities. Here, we show that the concept of contact can be generalized to an arbitrary partial- 
wave scattering. Contact of all partial-wave scatterings form a contact spectrum, which establishes 
universal thermodynamic relations with notable differences from those in the presence of s-wave 
scattering alone. Moreover, such a contact spectrum has an interesting connection with a special 
bipartite entanglement spectrum of atomic quantum Hall states, and enables an intrinsic probe of 
these highly correlated states using two-body short-ranged correlations. 


Ultracold atoms interact with short-range interactions 
U(r), which vanishes when the separation between two 
atoms is larger than a length scale ro- The atomic density 
in typical experiments is very low such that fcf < r 0 _1 is 
well satisfied. As originally discovered by Shina Tan (T|-[jl] , 
such diluteness leads to a wide range of fundamental re¬ 
lations among thermodynamic quantities in the presence 
of s-wave scattering. These universal relations are gov¬ 
erned by contact C. The momentum distribution n(k) 
has an asymptotic behavior C/|kj 4 at large |k|. The same 
C also shows up in the energy functional, the adiabatic 
relation, dynamic structure factors, and other relations 
BH3. Many of these universal relations have been ver¬ 
ified in experiments, and contact has been continuously 
inspiring physicists in both the ultracold atom and nu¬ 
clear physics community to explore its applications [8lil8| . 

In this Letter, we show that the concept of contact 
can be generalized to an arbitrary partial-wave scatter¬ 
ing, and we define a contact spectrum {Ci m } formed by 
contact in all partial-wave channels, where ( l,m) are the 
quantum numbers for the angular momentum. {Ci m } 
controls universal relations, such as the large momentum 
distribution and the energy functional, in the presence 
of arbitrary partial-wave scatterings, which have consid¬ 
erable differences from those in the presence of s-wave 
scattering alone, due to the fundamental differences be¬ 
tween high-partial-wave scatterings and the s-wave one. 
In addition to thermodynamic relations, { Ci m } provides 
physicists a new means to trace many-body physics from 
short-range quantities in strongly correlated systems. 

We consider the many-body wave function of a single¬ 
component atoms, either bosons or fermions, with a fixed 
total angular momentum L , 4' 3 £ ) = \H(ri, r 2 ,..., r,/v), 
where N is the total particle number. Due to the length 
scale separation ro <C kp 1 , T 3 £) has a unique asymptotic 
behavior when the distance between two particles r tJ is 
much smaller than kf, , 
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where r^ = r, - r j, R, ; = { r, + rj , r kjti,j}, E is the total 
energy, ipi m (r— rj ; e) is an unnormalized solution, which 
satisfies the boundary condition at r 7 = rj, for the rela¬ 
tive motion of the two particles with angular momentum 
quantum number (Z, m) and energy e under the two-body 
Hamiltonian, = — fjV 2 + Z7(|rj 3 -|). The symmetric 
or antisymmetric nature of the bosonic or fermionic wave 
function automatically picks up even or odd Z in the sum¬ 
mation. Equation Q can be understood from that other 
atoms essentially cannot affect the relative motion of the 
*th and jth atoms when |r,y| <C kp 1 . G; m (Ry; i? — e) 
characterizes the center of the mass of this pair of atoms 
and the other N — 2 ones, and acts as a “normalization 
factor” for the pair wave function — rUnlike 

an isolated two-body system, neither the angular momen¬ 
tum (Z,m) nor the energy e is conserved for any pair of 
particles in many-body systems. 

Whereas in the region r < r 0 depends on 

the microscopic details of the potential U(r), it takes a 
simple and universal form for r > ro where U (r) = 0, 


tpim('r l j-,e) ~ {ji(qerij) - teia[qi{q e )]ni{q e r io )}Ti m (r,y) 

( 2 ) 

where q 2 /M = e, ji(x) and ni(x) are spherical Bessel 
functions of the first and second kinds respectively, 
Yim (Tij ) is the spherical harmonics and ry 7 is a unit vec¬ 
tor. Using asymptotic forms of ji(x) and ni(x) at small 
x , we see that at small | r»j-1, 
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where co — _ UW _ R, — ( 2 Z— 2 s—i)i! r™ 
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phase shift has a standard low energy expansion, 

qf +1 cot[rn(q e )} = (-a z _1 + rfq 2 t + ...)rQ 21 , where a; and 

rf are the scattering length and the effective range of 

the Zth partial-wave [20]. For s-wave scattering in broad 

resonances, r§ is not important, and V , oo( r ij! e) has the 

conventional asymptotic form (y 1 ? - at small distance. 
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FIG. 1: (A-B), momentum distributions in unit of r(j of p- and d- wave two-body bound states respectively. Solid blue curves 
are exact results, dotted red ones represent the leading term fc 2! ~ 4 alone. Dash (pink) and dash dotted (green) curves include 
contributions up to the subleading and the third leading terms respectively. Insets compare the realistic momentum distribution 
n(k) and n(k) obtained from extending the universal behavior of n(k) in the regime k -C r,/ 1 to infinity in the zero-range 
interaction approximation limit. (C) {pi} for different bosonic QAH states of 4 atoms. With increasing the total angular 
momentum L, {pi } becomes broader and p;< m eventually vanish in the Laughlin state with filling factor 1/2. 


An angular-momentum-selective contact is defined for 
a given 1-th partial-wave as 


C lm = (4t t) 2 N(N - 1 )j dRij | J deGUR^E 

where f dR t j = f d( Ti 1 2 Tj )drk^ij- Using 


<0 , 

(4) 



dr^(r i,r 2 ,...r JV )e jk ' ri 


(5) 


a straightforward calculation as presented in the Supple¬ 
mentary Materials shows that 


tik = ^n im (fc)|F im (k)| 2 , (6) 

lm 

where the asymptotic form of ni m {k) is written as 

n lm (k) - 1 C lm k 2l - 4 + J2CLk 2l ~ A ~ 2s - (7) 

S —1 

In the presence of s-wave scattering alone, the expression 
reduces to the well known result rik —> Coofc _4 |Too(k)| 2 , 
where |F 0 o(k)| 2 = 1/(47 t) leads to a trivial pre-factor 
difference of Coo from the original s-wave contact de¬ 
fined by Tan. For any given high partial-wave scattering 
l > 0, the leading term in the large momentum distri¬ 
bution Ci m k 2l ~ 4 is distinct from the s-wave one. It also 
contains the contributions from the subleading and other 
terms, Cl m k 2l ~ 6 and etc. Expressions for Cf m are given 
in Supplementary Materials. Equation ([7]) can be demon¬ 
strated using a two-body problem, as shown in figure 
(1A-1B) for a simple square well interaction potential. 

Equation ([tJ) tells one that, if a zero range potential 
approximation r 0 —> 0 is applied, an ultraviolet diver¬ 
gence of energy shall occur. For instance, the lead¬ 
ing term k 2l ~ 4 in the Ztli partial-wave alone leads to 


f A dkk 2 ni m (k)Yi m (k) ~ A 2i+1 , where A is the large mo¬ 
mentum cutoff. For l = 0, Shina Tan has invented a 
remarkable functional integral form to remove such a di¬ 
vergence, which reveals the universal relation between 
the internal energy of dilute many-body systems and s- 
wave contact m- Despite that the scattering in a high- 
partial-wave has fundamental differences with the s-wave 
one, a similar energy functional could be derived in the 
presence of arbitrary partial-wave scatterings, 
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where E is the internal energy, n(k) is the momentum 
distribution by extending its universal behaviors in the 
regime k F k -C 1/ro to infinity, as shown in figure 

(1A-1B), and Ul = ^8^, and /f (kr 0 ) are 

a set of functions that depend on the range of interaction 
ro- The details of the derivation of the energy functional, 
as well as the explicit expressions of //(Zero), are given 
in the Supplementary Materials. In the regime Zero 1, 
for any s, //(fcr 0 ) —► 1. Using equations © and ([T]), one 
sees that the divergence in the kinetic energy encoded in 
n(k) is removed. On the left hand side of the equation 
a factor P ro =1 — N(N — l)/ ro /2 is included, where 



dciplmi^ij , ^)Glm(R'ij , E e)|~ 


(9) 

characterizes the total weight of the two-body wave func¬ 
tion in the volume defined by r,;,y | < ro- It can also be ex¬ 
pressed using Ci m (Supplementary Materials). We have 
also verified this energy functional using the two-body 
problem which can be solved exactly. 
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Whereas equation ([8]) is similar to the one for the s- 
wave scattering, it has a few notable features. First, 
the effective range rf explicitly enters the expression, 
since it is required for describing microscopic physics of 
a high partial-wave scattering, such as the bound state 
energy ~ l/(cqrf). Second, the terms Cf m must be in¬ 
cluded for completely removing the divergence for l > 0, 
where all terms in equation ([7]) contribute to the diver¬ 
gence. Finally, the singular part of the two-body wave 
function ~ r“/ _1 in equation |3| is not normalizable for 
l > 0 in the zero range potential approximation that 
takes r 0 — > 0. P ro is therefore required in equation ([8]). 
In the presence of s-wave scattering alone, the zero range 
approximation can be safely applied and P ro is negligible. 

All the above discussions can be directly generalized to 
two dimensions. One simply needs to replace ipi m (jij ; e) 
in equation (jTJ) by the solution of the two-body Hamil¬ 
tonian in two dimensions. For |zy| > ro, where Zj, = 
Xi + iyi is the coordinate of the 1th atom in two dimen¬ 
sions, such solutions are simply Bessel functions, i.e., ~ 
[cot i)i(q e )J(q e \zij\) - Yi(q e [zij\)\e tWi t, where ip(q e ) is the 
two-dimensional phase shift, and 6ij = arg {zij}. In the 
regime r 0 < \zij\ < kp \ 'F 2 d -t A 0 (Z y ) ln(^) +$^°, 



— 2Z + 1 

where A z (Zjj) = -T f 2t)! , JdeG^Z^-E - e) and 
Bi( Zij) = f(2 1- 2)UfdeG l (Z ij ;E- e), so that 

-2Z+1 

Ai(Zij) = Bi(Zij) ai 2 (2Z)!!(2Z - 2)!!’ ^ 

where Z tj = { , Zk¥=i j }. 

To simplify notations, we keep only the leading terms 
in the expansions of the Bessel functions. Similar to 
three dimensions, we define contact in an arbitrary high 
partial-wave channel, 

<Vo = pf-l)!!] 2 ^^ I dZ ^ B '( Zi ^ 2 ' ( 12 ) 

where f dZ.^ = f dZ l3 dzk- The leading term in 

the large momentum distribution is given by C)fc 2i_4 , 
which has the same power as that in three dimensions. 
Whereas the energy functional can be also derived in the 
same manner as that for three dimensions, we here focus 
on the application of contact spectrum in two dimensions 
in atomic quantum Hall states (QHS). 

QHS are intriguing many-body states carrying a large 
amount of angular momenta. At filling factor of 1/m, 
QHS is characterized by the Laughlin wave function 23 ] , 

= Y[(zi - Zj ) m e-^^ 21 ^ 2 /Vat, ( 13 ) 

i>j 


where A f is the normalization factor, m is odd(even) for 
fermions (bosons), a is the magnetic length, 'F [ J carries 
a total angular momentum L = mN(N — l)/2. There 
are two steps for establishing the applications of contact 
spectrum in atomic QHS. 

Step 1 We rewrite IF/, as 

m{2N-3) . |_. . | 2 \ 

Tl= £ p* FWZ^), (14) 

l—m ' 


where Zij = |Zi — Zj \e ieij , Xi = a I+1 V2 l+1 TTl\ so that the 
two-body wave function in the parentheses is normalized 
to 1, and Fs(Zij) is a set of orthogonal normalized wave 
function with angular momentum S and depends on (zt + 
Zj)/ 2 and the other N — 2 particles. 

Step 2 Though 'Fl is the exact ground state if one sets 
the interaction in all partial-wave channels with l > m to 
be zero in the framework of Haldane pseudopotential i!9j. 
it is in general a trial wave function for realistic systems. 
The exact ground state with the same total angular mo¬ 
mentum \Fff must contain corrections to *F£. For elec¬ 
trons with long-range Coulomb interaction, a standard 
approach is to numerically compute the overlap integral 
('FLl’iFff). In dilute systems, it requires that \Fff satis¬ 
fies the boundary condition at short distance in equation 
(10) and recovers \F l at large distance Zij \. 


Equation ([14| shows that each Laughlin state has a 
unique distribution of {pi}, with finite pi between l m i n = 
m and Z max = m(2N — 3). The very broad distribu¬ 
tion of {pi} as reflected by l m ax ~ lmin ~ N is a direct 
consequence of strong many-body entanglement in QHS. 
The same strategy can be applied to other types of QHS, 
such as the Pfaffian and quasi-hole states, described 
by = Pf ( 5 ^ 57 ) IW*- Zj)e~^\ 2 /^ and 


'Zh = 


a a 

dzi dz2 


d 

dz N 


0 -Y.M | 2 /4<r 2 


re- 


(n,^ - e 

spectively, where Pf is Pfaffian of (zj — Zj)~ l ■ 

Figure(lC) shows {pi} for a few bosonic QAH states. 

Equation (14) also allows one to define a special bi¬ 
partite entanglement spectrum. For a pure state wave 
function, |\F) = ^ where Z = 

Yf e~^ x , and |\Fj*) and |\F^) are the eigenstates of its two 
subsystems, {£a} defines the entanglement spectrum [22]. 
Using equation ( [l4| , one observes that {£; = —21ogp/ — 
logZ} characterizes the entanglement spectrum between 
the relative motion of a pair of atoms and the rest of 
the system in the angular momentum space. Whereas 
entanglement is usually considered by spatially dividing 
the system into two parts, {p{\ here provides one a new 
angle to characterize QHS, since {pi} captures the corre¬ 
lation in the angular momentum space, regardless of the 
distance \zij\ between particles. Thus, the signature of 
many-body entanglement persists in short distance and 
is captured by the contact spectrum as discussed below. 
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When \zij\ approaches zero, equation (14) becomes 


rn{2N-3) 

-> 51 K l Pi F L-i(^j)\zij\ l e iW ' J , ( 15 ) 

l—m 

which does not satisfy the boundary condition specified 


by equation (10). Though the high partial-wave scatter¬ 


ing length in general is smaller than the s-wave one, for 
realistic interaction with a finite range r 0 , the singular 
terms l/|z»j|* must exist and becomes important in short 
distance. In ultracold atoms, aj^o can be tuned by Fesh- 
bach resonance and other techniques, so that the singular 
terms could have observable effects even away from the 
resonance. 

Whereas i,}®*^) can also be computed for short- 
range interactions, the diluteness provides one a simple 
method to qualitatively estimate such an overlap inte¬ 
gral. We define a length scale r* = max{r ; *}, where 
r* satisfies \Ai\rf l = \Bi\r}~ 1 . It is easy to verify that 
r* ~ ro(ai/ro) 2 i for low energy scattering « rj 1 . If 
the separation between any paired atoms is much larger 
than r*, the correction to the wave function due to the 
weak scattering cii <C cr in all partial-wave channels be¬ 


comes negligible. Using equation (10), the exact wave 
function becomes 


qj f ‘ r x 


r*«|ay|«<r 

-> £>( Z «) | 

Iji 0 


za\ l e iWi - 


(16) 


Compare equation (15 16), one see that, if the follow¬ 
ing equation 


A i(Zij) = A ; 1 piF L ^i(Z i j), 


(17) 


is satisfied, the exact ground state wave function 4'“' re- 

| Zij |]^>r* 

duces to 'S’l at large inter-particle distance, 4>f / x —► 

\Fi. The overlap between 4/f x and 4 >l can be estimated 
as 1 — 0{k 2 F r* 2 ), where ~ 1 in QHS. The criterion 
for 4/^ to be a good approximation of the exact ground 
state is that r* <C kp 1 . If the scattering for all l > m van¬ 
ish, rf >m = 0, and 4becomes the exact ground state, 
consistent with results from Haldane pseudopotential. 

Using equations ( 11|12|17 ), one establishes an intrinsic 
relation between contact spectrum {Ci} and the entan¬ 
glement spectrum {p/}, 


Pi = 


l X i 


2 l+2 H ^N(N - 1) 


c; 


1/2 


(18) 


As contact spectrum can be measured through various 


schemes, equation (18) allows one to probe {pi}. Besides 
the large momentum distribution measured in Time-Of- 
Flight experiments, another useful scheme for measuring 
{Ci} is photoassociation, which can be made angular- 
momentum-selective. For instance, by tuning the laser 
to be resonant with an excited state with a particular 


angular momentum l e ^ 0, only the l = l e th partial- 
wave of the relative motion of two atoms contributes to 
the photoassociation. For a two-body problem, an im¬ 
portant quantity to characterize the stimulated rate of 
populating the excited state and the scattering lengths 
of the optical Feshbach resonance is the coupling strength 
to the excited state, r[ 2 ^ = c~ 1 (2nl)d 2 , 1 /|d>j 2 h| 2 [23V 

S3] , where c is the speed of light, / is laser intensity, 
is the dipole moment, ip e y is the electronically excited 
molecular wave function with angular momentum Z, and 
(f>\ is the relative wave function of two atoms. 

In dilute many-body systems, the size of the electroni¬ 
cally excited molecules is much smaller than inter particle 
spacing. Since G;(Z,j; E — e) (or — e) in three 

dimensions) acts as the normalization factor for the two- 
body wave function the coupling strength in 

the many-body system can be written as 

| 2 , (19) 


where a low-energy expansion ipi(zij,e) = X) s £S V’f ( z b) 
has been used. In an experiment by Gemelke, et al, pho¬ 
toassociation has been used for detecting atomic QHS 
clusters [26]. Whereas the decreased signal of photoasso¬ 
ciation with increasing the rotation frequency is consis¬ 
tent with accessing the QAH regimes, an angular momen¬ 
tum selective photoassociation can further probe the rise 
of high partial-wave contact, so that much richer infor¬ 
mation of atomic QAH states can be obtained, as shown 
by figure (1C) and equation (18). 

In conclusion, we have generalized the concept of con¬ 
tact to arbitrary partial-waves, and defined contact spec¬ 
trum {Cim} and {Ci} in 3D and 2D. In addition to ther¬ 
modynamic relations, we have applied {Ci} in atomic 
QHS, and shown that it serves as an intrinsic probe of 
atomic QHS due to a connection with a special entan¬ 
glement spectrum of QHS. We hope that this work may 
stimulate more studies on the connection between short- 
range correlations and many-body physics in strongly 
correlated dilute systems. 

Note: near the completion of this manuscript, two the¬ 
oretical and one experimental paper ¥I7\W3\ addressing 
contact of a single high partial-wave scattering (Z = 1 
case in our work) showed on arxiv. All works agree 
on the leading term in the large momentum distribu¬ 
tion for the p-wave scattering, and [28] [29] have also 
explored the subleading term consistent with ours. QZ 
acknowledges D. Wang, J. H. Thywissen, S. Zhang and 
J. Zhang for discussions. This work is supported by 
RGC/GRF(14306714). 
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Supplementary Material 


In this supplementary material, we present the results on the large momentum distribution and the energy func¬ 
tional. 


Large momentum distribution 

.2 

Starting from n k = V]. j 11 ,,-., dr, / drj\i>(ri, r 2 , ...rjv)e lk ri , we define ^(k) = / dr,\I>(ri, r 2 , ...rjv)e* k ' ri so that 
n k = f U jf H ^ r j IAi( k )| 2 • We also define a purely two-body quantity, 


Am(k;e) = rj) ^ 


( 20 ) 


so that in the regime |k| » kp can be written as ^(k) —> Y2j& e lk ' r;; Y2im f deG; m (Rjj; E — e)Fi m ( k;e). If one 
extends the wave function in the regime [ro,oo] to [0, oo], ,F(k;e) in the regime kp <C k is given by, 


Fim{ k; e) 


For s-wave scattering, we have 



Al9e 

l-l 



F 0 o(k;e) = 4 t x/k 2 


For l > 0, using e lk r = 47 t]C^ 0 Yl m =-i * Z A(fcr)Y im (k/fe)Y'^(r/r), we obtain 


( 21 ) 


( 22 ) 


V>i( k ) 


>2 < ,k 


Z^ 


47 xi' 


( k l ~ 2 


deGi m ((R.ij; E — e) + k / deq e Gi m (R,ij; E — e) + ... )Y; TO (k) 


(23) 


When calculating the integral / ^ r jl^i( k )| 2 ! cross term e * k ( r j' — r j) vanishes in the large k limit. Changing 

dvj^ti to dUij, the cross term G), m , (R,; 7 -, E — e q )Gi m (Rij, E — e q >) also vanishes, due to the orthogonality of wave 
functions with different angular momenta, we obtain 

n k = fc -4 Cb 0 |loo(k)| 2 + £ {k 2l - A C lm + k 2l ~ Q Cl m + k 2l ~ 8 C 2 m + ...)|Y Zm (k)| 2 , 

1^0,m 


( 24 ) 
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Cim has been defined in the main text, and 

Cl m = 32 tt 2 N(N - l)Re\ / dR„ ( / deG^Ry; E - e) 


J deq^Gim (R, j ; E — e 


(25) 


Qm = 32n 2 N(N — 1) J dRy{ite[^y deGJ^Ry; E - e)) (^J deq 4 e G lrn (R lj; E - e)j ] | J deq 2 e G lm (R i:j] £-e)| 2 }. 

(26) 

Other C* m can also be written down straightforwardly using the same procedures. 

Energy functional 

In the regimes where the distance between any pair of particles is larger than r o, the many-body Schrodinger 
equation becomes ^ = where the arguments (ri, r 2 ,rjy) has been suppressed. It leads to 


n«u- S '’’l*= s / n**i*i 


(27) 


where ]+ dr,; means the integration is carried out in the regions where any pair of particles is larger than r o. 
In the zero range interaction limit, one extends the asymptotic form 




J de ( tan q [m(q e )f l ^ 9eri ^ ~ 9 ' +ln, ^ er<J ^) G i™(Rij-,e)Yim(rij), 


(28) 


which is valid in the regime r o <C r. (J < to the regime < ro- Define <f>, which is identical to if 't'ij > ro 
and given by the right hand side (RHS) of the above equation, the integral on the left hand side (LHS) of equation 
(27) can be extended to the whole real space, and meanwhile, the integral in the unphysical region r,j < ro must be 


subtracted to remove the divergence. If three-body physics is ignored, we obtain, 

t2 ^ r tt r 

2 M 


lhs=i n *<*• i £ w - / n ^ r ‘ +r< 


r „ h 2 k 2 _ iv(iv -1) 

= / dk _ n k — 


2 


2 M 


dR,: 






/o 


I $ 


2 ./ 


f „ h 2 k 2 N(N - 1) Ti 
= I dk^^n k + v - ’ 


V 2 - 

V tj 


2 M 




M 


^ dr fj{[ F Zm(%) J deGi m (R tj ;e) ( tan ^ g ) ~ g+nz+r^) 


rv^ + 1 


tan[r7;(g c )] 


jliQenj) - q[ +1 ni(q t r i: j) j Yi m (r tj ) | 


(29) 


As seen from the asymptotic form of n(k) at large |k|, we conclude that the integral f dk r ^jj-fik is divergent. It 
is worth mentioning that not only the leading term fc 2 ^ 4 contributes to the divergence, many other terms for large l 
will also contributes to the divergence. The divergence therefore needs to be carefully removed. 


One term in the last line of equation (29) automatically removes the divergence. To see this fact, we examine the 


contributions separately from ji{q e rij) and nz+ry,-). Since V+f+r+Fz+rj,,) = -qiji{q e .r%j)Yi rn {T i j), one sees that 

7+1 


dr, 


/{[*£,(*«)( 


qi -\-1 

jMeVij) - q[ +1 m{q e rij) 


\tan[? 7 /(g e )] 


Vy tan [rn(qe)\ JlyHt ' lJ 


jl (QeT'ij)Yim (** ij ) 


2J.+1 

{<lt.ro) 


tan[ry ; (g £ )] 


(30) 

In the low energy limit, (q e ro) 2 •C 1, such a contribution is negligible. 

In contrast, ^ 2 j ni{q t r i j)Yi m {i i:j ) is crucial. By making use of ^£r 2 £ - +r = -(2 1 + 1 )S(r)/r l+2 

and 


2Z+1 


tan [rn(q )] = ( — + r iQe + • ••) r o 2 > an d expansion of ni(q e rij) at small one could compute 

fo° ^ / im(^j) n l(^ r iy)^ij n l(9e r ij)Yim(rij) systematically. Not surprisingly, this removes the divergence in J dk^A-h k , 
since the divergence is indeed caused by the singular behavior of rcz(g<Ty) at small distance. 
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As a demonstration, here we show how to treat / = 1, where we have 


N(N - 1) h 2 
2 M 

N(N-l) n 2 
2 M 


E / dR o- r°rfr iJ {[y 1 * ro (f i ,)E 


/3i s ffr 


,*2s 


2—2s 
s =0 U 


[s 


Pi»9l 


2 s 


1 V 


2-2s Ilrn 

's =0 


E / dR *. 


dr, 






r 2 

u 


+ Pn9u 


Pl09lm 


(*«)]} 

Pll9lm 


Y, 


lm 


(*«)} 


and 9im = I deq 2 e s Gi m ( R tJ ;e). 
The leading term gives 


zv(.zv — i) h 2 

2 M 


E / dR t 


dr 




_ ^ r dk a 2 * 2 j- c 

(2tt) 3 TM l 


lm 


(2tt) 3 2M l k 2 


K, 


lm 


/ W* ^ \ u2 /^lO^lm 

(k)y lm (k)| 


*J „2 

ij 


The subleading term gives 

N(N — 1) ?i 2 
2 M 


E / dR ^ 


dr^-y, 


U J lm 


(*«){ 


Piog\° 


V 2 ,/3n5? m + ^n5imV 


fiw 9° 


lm 


iJ r 2 
ij 


J* ^lm (^r?) 


=E 


dk ^ 2 fc 2 r Cj. 


(2tt) 3 2Af 


{ - -jrfi(kro)Y? m (l)Y lm (l)} 


where fl(kr 0 ) = 1 - |jo(fcr 0 ). 
Similarly for l = 2, one has 


A'N = i - 


fi(hn) = 1 - + jo(fcro)] - 


> ji(fcrp) 
kr a 


Ct - (47r) 2 7V(7V - 1) / dRy^™ 

For any Z, there remains a term Y^ n (r ij )jf(q € r i j)V^ j ni(q t r i j)Yi rn (r ij ). 
Again, by applying = _f97 xiUW/,1+2 


2 dr 1 dr 
2 


-j = —(21 + 1 )5(r)/r l+2 , we see that such a term leads to 


N(N- 1) ft 
2 M 


E / dR^^d^lfy^lf^) J deGU Ry;e) 


' \tan[» 7 i( 5 r e )]- / V J ‘ 

deG/ m (Rij, e) ( q € l fii(q e Tij')') j* 


= E 


t-2 -27 

n r n 


C)m 


a 2 - 22 


32tt 2 M a; 64 tt 2 M 


0 e^l 
V-;, 


n ^lm 


(31) 


(32) 


(33) 


(34) 

(35) 

(36) 


(37) 


1 # Q ~ 

The reason that only Ci m and Cf m show up here is because the low energy expansion of the phase shift ta E( g y = 
(—A + rfg 2 + ...)7 -q 21 has been used. 

For the right hand side (RHS) of the equation (27), since is normalized to 1 and three-body physics leads to high 
order contributions in the dilute limit, we have 


ndr^| 2 = l- 


N(N — 1) 


E J dR b J dr*, J deipi m i+ij i e)G l m (R ij , E — e) | z 


In low energy limit, we can expand the pair wave function e) as: 

1>lm(rij-,e) -t E^ S) ( r «)9? a 


(38) 


( 39 ) 






































The equation (38) can be expressed by using contact: 


Pr o = l- 


N(N — 1) 


E/r° dr bl f *(EC S) (^)? £ 2s )Gz m (R U) i?-e)| 2 

Zm ^ s 


= i-E 




L32tt 2 


lw™( r «)l 2 *ij + ^ 


'lm * ( r n ) Wm ( r « ) + Wm * ( r « Mm ( r t? )) dr 


(40) 


All the integrals are purely two-body quantities, which can be obtained by either microscopic calculations or 
measurements in simple many-body systems where C/ m are known. 







